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Abstract. Rare B meson decays are an excellent probe for beyond the Standard Model physics. Two very 
sensitive processes are the b — > and B° sd — > decays. We report recent results at a center of mass 

energy of -\fs = 1.96 TeV from CDF II using 7 fb -1 at the Fermilab Tevatron Collider. 



1 Introduction: b -» decays 

b — > sp + p~ decays are flavor changing neutral current 
(FCNC) processes that can only occur through higher or- 
der box or penguin amplitudes in the Standard Model. New 
physics can be probed by measuring various combinations 
of their decay rates. One of the most sensitive observables 
is the forward-backward asymmetry of the muons (Afb) as 
a function of the squared di-muon mass. 

Decays of interest are B + p + p~K + , B° d ^ p K*°(892)., 
B° s n + M -c/>, B + p + p~K*\ B° d p + p~K s , A,, -» 
p + p~ A. The CDF analysis that uses 6.8 fb of data is de- 
scribed here Q, Q- 



1.1 Analysis Method 

1.1.1 Branching Ratio Measurement 

Branching ratios for the b — > sp + p~ decays are measured 
relative to normalization modes, where the two muons orig- 
inate from a J/tf/ decay. For the event reconstruction CDF 
requires two muons with a transverse momentum (pj) greater 
than either 1.5 GeV/c, 2.0 GeV/c, or 3.0 GeV/c depending 
on the muon trigger. The six modes are then reconstructed 
where the K*°(892) is reconstructed from #*°(892) -> K + n-, 
the (p from K + K~, the K* + from K* + K s n + , the K s 
from K s — > n + n~, and the A from A — > pn~ . To avoid con- 
tamination from resonant modes such as the J I if/ and ij/' , 
candidates with di-muon masses near these resonances are 
rejected. 

The events then have to meet loose preselection re- 
quirements before an artificial neural network (NN), which 
combines multiple discriminating variables into one vari- 
able, is applied. Signal is modeled with p^-reweighted Pythia 
signal Monte Carlo simulations (MC). The reweighing is 
done by comparing the MC pj distribution with that of the 
normalization modes. The background is modeled by sam- 
pling the b hadron (Hi,) mass sideband regions. 

The final signal yield is obtained by an unbinned max- 
imum likelihood fit to the Hi, mass distribution. The proba- 
bility distribution function (PDF) of the signal is parametrized 



with two Gaussian with different means while the the back- 
ground PDF is described by a first or second order poly- 
nomial. Peaking background contributions are subtracted 
from the fit results for the signal yields. The only signifi- 
cant peaking contribution is cross-talk among B'j — > p + p~K* 
and B° s — > p + p~<p which has a ~ 1 % contribution to the total 
observed signal MC yields. The final branching ratios are 
calculated as follows: 



<B(H h -> hp + p~) 
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where Hb signifies the Z?'j, B + , B° s , or Ai, and h represents 
#*°(892), K + , K* + , K s , A, or 



yields after the optimal NN selection, 



and M°° s f are the 

.loose 

is the relative 



efficiency of the loose selection cuts, and e^J _ is the NN 
cut efficiency on the loosely-selected events. The NN is 
not applied to the normalization mode because the sig- 
nal/purity and size is sufficient with the loose selection 
cuts. The NN cut efficiency are obtained from signal MC. 

The three leading systematics are the systematics on 
the efficiency, S(Hb — > Jjiili), and background PDF. The 
main sources within the efficiency systematic errors are the 
MC reweighing and the the trigger turn on. The total sys- 
tematics for the B + p + p~K + , B° d -> p + p 7T°(892), 
B° ? p + p-(f>, B + p + p-R* + , B° d -> p + p-R s , A b 
p + p-A are 5%, 6%, 32%, 8%, 6%, and 32%, respectively. 
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1.1.2 Forward-Backward Asymmetry Measurement 

For the B d p + p-K*°(%92) and B -> p + p~K* (B d -> 
p + p-K*°(%92) and B + -> p + p-R* + combined fit for in- 
creased sensitivity) decays the forward-backward asym- 
metry (AfB), the K* longitudinal polarization (fY), the trans- 
verse polarization asymmetry (A^ ), as well as the triple 
product asymmetry of the transverse polarizations (A Im ) 
are measured. For B + — > p + p~K + only Aps is measured. 
The A fb and Fi angular quantities are extracted from cos 6^, 
the cosine of the helicity angle between the p + (p~) mo- 
mentum vector and the opposite of B (B) meson momen- 
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turn vector in the di-muon rest frame, and cos 9k, the co- 
sine of the angle between kaon momentum and the oppo- 
site of the B meson momentum vector in the /T"°(892) rest 
frame. and A Im are extracted from <p, where the angle 
<p is the angle between the two decay planes of the di-muon 
pair and K — n decay pair. B decay amplitudes are calcu- 
lated using operator product expansion and Wilson coeffi- 
cients. There are many non-SM predictions for Af B from 
models with different Wilson coefficients 0. 

The angular measurements are extracted using an un- 
binned maximum likelihood fit containing B mass shape 
signal and background PDF's as well as signal and back- 
ground angular shape PDF's. The mass shape PDF's are 
divided into several di-muon mass bins and are described 
as 25 bin histograms. The combinatorial background PDF 
is taken from the B meson higher sideband. The angular 
acceptances are also described as 25 bin histograms and 
are derived from phase space signal MC. 

As a control A FB and F L are fitted to B° d -> J/t//K*°(S92) 
and AfB only to B + — > J/iffK + . This cross check yielded 
measurements that were consistent with other measurements. 

The dominant systematic uncertainties for all angular 
measurements are the signal fraction and B mass shape 
uncertainties. This uncertainty is assessed by varying the 
signal fraction and shape parameters by + lcr (statistical 
uncertainty) when performing the B mass fit. 



1.2 Results 

The resulting yield for the six decays are shown in Fig- 
ure Q] All measured branching ratios agree with previous 
measured values as well as theoretical predictions. CDF 
reports the first observation of A/, — > pfpC A with a signifi- 
cance of ~ 6<x and a measured branching ratio of S(Ai, — » 
H+^-A) = (1.73 + 0.42 [staf] +0.55 [syst]) x lO" 6 . 

The results of the angular measurements for the com- 
bined B — > n + n~K* fit are shown in Figure [2] They are 
compatible with theory predictions and competitive with 
the results from the B-factories. 
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2 Introduction 



s,d 



fi + fi" decays 



Fig. 1. Yields for the various b — > sn + [i decays. 



2.1 Analysis Method 

CDF collects opposite sign muon candidate events using 
di-muon triggers. 

B + — » J/if/K + events are collected on the same trig- 
gers as a relative normalization mode to estimate S{B° S — > 
yU + //~) as: 



B® d — > ju + yu are FCNC decays that are highly suppressed 
by the SM. The SM predictions for these branching frac- 
tions are S(B° S -> fx^/T) = (3-2 ± 0.2) x 10~ 9 and S(B° d -> 
// + yu") = (1.00 + 0.1)x 10~ 10 [!4]. These predictions are one 
order of magnitude smaller than the current experimental 
sensitivity. Previous bounds from the CDF collaboration, 
based on 3.7 ftT 1 of integrated luminosity, are S(B° S — » 
fffT) < 4.3 x 10~ 8 and <B{B d -> n + fi~) < 7.6 x 10 l9 at 
95%C.L. 0. 

Enhancements to the rate of B° s — > decays occur 
in a variety of different new-physics models. In supersym- 
metric (SUSY) models, new particles can increase S(B° S — > 
ju + yu _ ) by several orders of magnitude at large tan/?. In the 
minimal supersymmetric standard model (MSSM), the en- 
hancement is proportional to tan 6 /?. Global analysis includ- 
ing all existing experimental constraints suggest that the 
large tan/3 region is of interest |6 7I8I9L 

This document describes the current status from CDF. 



S(B° S -» jiV) = 
■S(B + 



J/t/fK + ) ■ 8(J/if> 
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where Ngo is the number of candidate B ( ' — > events, 
a B o is the geometric and kinematic acceptance of the di- 
muon trigger for B® — > ^i + ^r decays, e^ ml is the total ef- 
ficiency (including trigger, reconstruction and analysis re- 
quirements) for B° s — > y^yT events in the acceptance, with 
Nb+, a B +, and e™ al similarly defined for B + -> J/tj/K + 
decays. The ratio f u /f s accounts for the different fo-quark 
fragmentation probabilities and is (0.402 + 0.01 3)/(0.1 12 + 
0.013) = 3.589 + 0.374, including the (anti-)correlation be- 
tween the uncertainties |10|. The final two terms are the 
relevant branching ratios T>{B + — > J/iffK + ) ■ &(J/ift — > 
//V) = (1.01 + 0.03) x 10~ 3 ■ (5.93 + 0.06) x 10~ 2 = 
(6.01 + 0.21) x 10~ 5 CI. 
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Fig. 2. Results of angular measurements for B — > K*. 



The CDF analysis described is also sensitive to B° d — » 
decays. The value of &{B° d — > is estimated 

from Equation[2]substituting Z?'j for Z? ( ', and changing /„//? 
to f u /fd = 1. All other aspects are the same as the B° s — » 
search except where noted below. 

The analysis is done by first estimating the acceptances 
and efficiencies, then creating a multivariate discriminant 
for background rejection. This discriminant is optimized 
a priori with Pythia signal MC and data mass sideband 
events and validated with the B + — > J/i(/K + normalization 
mode. Before the discriminant is applied the events are re- 
quired to pass baseline cuts that consist of loose require- 
ments on muon ID, vertexing related variables, muon pj, 
and di-muon pj. The background is then estimated, which 
has two sources: combinatorial background and peaking 
background (B— > h + hr). Finally, when the background is 
well understood the di-muon mass signal region is unblinded. 



2.2 Signal and Background Properties 

The signal candidates are fully reconstructed events with a 
secondary vertex due to the long lifetime (~ 450pm) of the 
B { ] meson. Signal events will have a primary-to-secondary 
vertex vector that is aligned with the B° s candidate momen- 
tum vector. Another property of signal events that is unique 
is that they are very isolated (with few tracks near the muon 
tracks) due to the hard B fragmentation. 

Background events tend to be partially reconstructed 
and be shorter lived than signal. They also have a softer pj 
spectrum, higher activity of tracks, and misaligned primary- 
to-secondary vertex and momentum vectors. The combi- 
natorial background consists of sequential semi-leptonic 
decay (b — » cyTX — » yfyTX), double semi-leptonic de- 
cay (bb — » pTp + X), continuum di-muon events, as well 
as fake+p and fake+fake events. The expected number of 
combinatorial background events in the signal window is 
estimated by extrapolating the number of events in the side- 
band regions to the signal window using a first order poly- 
nomial fit. 



The peaking background from two-body hadronic B 
decays is also evaluated. These backgrounds are about a 
factor ten smaller in the signal region than the combi- 
natorial background and need to be estimated separately. 
The dominant contributions to this source of background 
are the decays of B° s and B° to h + h'~ final states, where h 
or h' can be either 7r ± or A' ± and are misreconstructed (fake 
muons). 



2.3 Analysis Improvements and Signal 
Discrimination 



CDF has updated their search with an improved analysis 
and significantly more (+3.2 ftT 1 ) data. The forward muon 
acceptance (left Figure has been increased and the neu- 
ral network has been improved to achieve twice the back- 
ground rejection (right Figure |3). Additionally the neural 
network has been extensively tested for mass bias. The 
peaking background prediction has also been improved. 
The central-forward (CF) channel of the analysis has in- 
creased in statistics by ~15% resulting in a total increase 
for both channels of ~7%. 
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Fig. 3. Left: CDF Muon detector layout. The blue region has 
been added for the updated analysis. Right: Comparison of signal 
and background efficiency for new and old NN. 



The final NN network consists of 14 input variables. 
Extensive cross checks were done to demonstrate the 14- 
variable network does not sculpt the di-muon invariant mass 
distribution. Figure|4]shows the correlation across the M^^- 
mass range. No significant correlations are seen. 
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Fig. 4. Left: NN output as a function of di-muon mass for the 
signal sample. Right: NN output as a function of di-muon mass 
for the OS- control sample (see Sect. l2.4b . 
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2.4 Control Samples for Background Estimates 

The background estimate procedure is cross checked on 
four independent background samples: 

OS-: opposite-sign muon pairs, which pass the baseline 
and vertex cuts and have a negative lifetime; 

SS+: same-sign muon pairs, which pass looser baseline 
and vertex cuts and have a positive lifetime; 

SS-: same-sign muon pairs, which pass looser baseline and 
vertex cuts and have a negative lifetime; 

FM+-: opposite-sign fake-muon pairs, at least one leg of 
which is required to fail the muon ID requirement, pass- 
ing looser baseline and vertex cuts for both positive and 
negative lifetimes. 

These are representative of various background con- 
tributions and thus are excellent control regions to test the 
method of background estimation. The estimated backgrounds 
are compared with actual observation. The agreement be- 
tween predicted and observed number of background events 
in the signal region is good across all control regions. 



2.5 Results 



The observed events for the B°. 



search are shown 



in the top part of Figure \5\ The data is consistent with 
the background prediction and yields an observed limit of 
S(B d -> /z + /T) < 6.0 (5.0) x 10~ 9 at 95% (90%) C.L. 

An ensemble of background-only pseudo-experiments 
are employed to estimate the significance as a p-value. 
The effects of systematic uncertainties are included in the 
pseudo-experiments by allowing them to float within Gaus- 
sian constraints. The resulting background-only /9-value 
for the — > /i + /i~search is 23.5%. 

In the B® search region the data exceed the background 
prediction in bins with AW > 0.97. The bottom part of 
Figure [5] contains the detailed breakdown of the expected 
background and actual observation for the individual NN 
and mass bins for the Z?" — > y^pT search. 

The source of the data excess in the 0.970 < AW < 
0.987 bin of the B Q S signal region was investigated. The 
same events, same fits, and same methodologies are used 
for both the B® and B°, searches. Because the data in the B°. 

s d d 

search region shows no excess, problems with the back- 
ground estimates are ruled out. Problems with the NN are 
ruled out by the many studies performed. The most plausi- 
ble remaining explanation is that this is a statistical fluctu- 
ation. 

The /7-value for the background-only hypothesis is 0.27%. 
Additionally a p-value of 1.92% is also calculated assum- 
ing the SM+background hypothesis. 

A Ay 1 fit is used to determine the S(B° S — > most 
consistent with the data in the B ( ' search region. A central 
value and 68% C.L. is calculated at S(B° -> /i + /T)=(l .8^)X 
10~ 8 . Additionally a bound at 90% (95%) C.L. on the branch- 
ing fraction of B° -> ju + >l/ is set at 5.4 x 10~ 9 < S(B° -> 
H + H~) < 3.9 x 10~ 8 (2.8 x 10~ 9 < <B(B° S -» n + n~) < 
4.4 x 10~ 8 ). 

Finally, upper limits at 95% (90%) C.L. of S(B° S -> 
fi + fi~) < 4.0 x 10~ 8 (3.5 x 10~ 8 ) are set with the CLs 
methodology. 
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Fig. 5. Top: Background estimates (light gray), systematic er- 
ror on background estimates (hashed area), Poisson error on the 
mean (error bars on points), and data for the B° d signal window. 
Bottom: Similar plots for S° with the addition of the SM expec- 
tations (dark gray). The data is divided into 8 NN bins, of which 
lowest 5 NN bins are combined into one bin for both figure, 5 
mass bins, and two muon topologies (CC and CF). 



3 Conclusion 

FCNC decays are a powerful probe for New Physics and 
CDF continues to lead the searches in the B sector. The 



analysis has yielded the first two sided bound 
H + jJT) and will be updated with the full CDF 



to %(B° S -» 

Run II data set. CDF also has the first observation of A\, — > 
H^pT A as well as angular measurements that compete and 
agrees with the B -factory results. A® and A Im are also 
measured for the first time for B — > p^yTK* decays. 
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